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I.  INTRODUCTION  -  SUMMARY  OF  THIRD  QUARTER  PROGRESS 

^  ■ 

This  report  summarizes  the  work  accomplished  during  the  third  quarter 
of  the  Spacecraft  Electric  Generating  and  Propulsion  System  Integration  Study 
Program  under  Air  Force  Contract  No.  AF33(657)-8488 .  The  primary  objectives 
of  this  program  are: 

1.  To  develop,  demonstrate,  and  provide  a  working  manual  and 
digital  computer  program  for  the  system  design  analysis  and 
optimization  procedi^re  known  as  LEADER. 

2 .  To  apply  this  procedure  to  space  power-propulsion  systems  em¬ 
ploying  the  SPUR  350  KW  and  1  MW  nuclear  space  power  plants. 

The  initial  phase  of  the  contract  effort  has  been  devoted  to  the  solution 
of  two  Illustrative  problems,  (1)  optimization  of  a  nuclear  space  power  plant- 
transmission  line  system,  and  (2)  optimization  of  a  space  power  plant  condensing 
radiator  p)anel,  as  a  means  for  satisfying  the  first  objective.  The  previous 
quarterly  prpgress  report  contained  the  results  of  the  model  formulation  analyses 
of  the  two  Illustrative  problems  and  the  results  of  the  optimization  of  the  trans¬ 
mission  line  problem.  The  optimization  results  were  obtained,  however,  with 
a  transmission  line  model  which  lacked  certain  of  the  parameter  interaction  terms. 
This  has  been  rectified  during  the  third  reporting  period  and  revised  optimization 
results  have  been  obtained.  The  radiator  panel  optimization  has,  in  addition, 
been  completed.  The  results  obtained  for  the  two  illustrative  problems  contain 
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(1)  the  values  of  each  of  the  independent  design  parameters  which  minimize 
the  over-all  system  weight  within  the  boundaries  and  constraints  specified, 

(2)  the  value  of  the  minimum  over-all  system  weight,  (3)  the  boundaries  or 
constraints  which  prevent  further  reduction  in  over-all  system  weight ,  and 
(4)  the  sensitivity  of  the  over-all  system  weight  to  small  changes  in  the 
optimum  design  parameters. 

The  process  of  solving  the  two  illustrative  problems  has  served  to 
clarify  the  capabilities  of  the  two  computer  programs  comprising  LEADER  and 
to  identify  the  various  procedural  elements  which  must  be  performed  in  order 
to  achieve  effective  utilization  of  the  LEADER  technique.  These  procedural 
elements  are  being  described  in  a  User's  Manual  in  order  to  facilitiate  the 
application  of  LEADER  to  subsequent  problems.  During  the  subject  reporting 
period,  the  first  part  of  this  User's  Mahual  has  been  completed  and  is,  therefore, 
included  as  a  supplement  to  this  progress  report.  This  section  contains  the 
procedures  recommended  for  the  model  formulation  phase  of  LEADER  problems 
and  Includes  wdiatever  tabular  data  may  be  needed.  The  results  of  the  model 
formulation  phase  will  be  a  series  of  analytical  models  of  each  of  the  sub¬ 
systems  comprising  the  over-all  system  under  investigation.  These  models 
will  then  be  assembled  into  an  over-all'  system  model  and  processed  by  the  optimi¬ 
zation  program.  The  second  section  of  the  User's  Manual  will  contain  the  procedures 
necessary  for  this  final  phase  of  the  analysis.  This  section  has  been  outlined  and 
detailed  definition  is  in  progress.  It  is  expected  that  this  Sfi,cond  section  of  the 
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manual  will  be  completed  within  the  next  few  weeks. 

The  process  of  generating  and  accumulating  basic  design  data  on  the 
major  sub-systems  for  the  350  KW  arc  jet  and  the  1  MW  ion  engine  power- 
propulsion  systems  has  been  continued  during  this  reporting  period.  Design 
procedures  and  performance  characteristics  have  been  established  for  the 
turbine,  and  the  generator  sub-systems.  These  procedures  will  be  described 
in  more  detail  in  the  following  sections.  The  task  of  formulating  the  analytical 
model  of  these  sub-systems  has  been  initiated.  It  is  expected  that  the  entire 
model  formulation  Process  will  be  completed  by  the  second  week  of  February. 

IL  ILLUSTRATIVE  PROBLEMS 

The  two  illustrative  problems  have  been  structured  to  serve  as  a  vehicle 

for: 

(a)  developing  the  procedural  elements  of  the  LEADER  technique  which 
would  be  required  to  supplement  the  existing  computer  programs, 

(b)  as  a  demonstration  of  the  type  of  information  required  for  input  and 
the  type  of  results  to  be  expected  from  the  use  of  LEADER  on  complex  system 
design  problems ,  and 

(c)  as  a  guide  in  the  development  of  a  User's  Manual  v^^ich  will  facili¬ 
tate  the  use  of  LEADER  for  such  problems. 
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The  mechanics  of  the  model  formulation  process  were  discussed  in 
relation  to  the  two  illustrative  problems  in  the  previous  progress  report  and 
will,  therefore,  be  omitted  from  this  discussion.  This  section  will,  instead, 
report  on  the  process  evolved  for  feeding  the  system  model  into  the  optimization 
program  and  on  the  type  of  results  obtained.  It  should  be  noted  that  the  number 
of  Independent  variables  and  non-linear  constraints  required  for  each  of  the 
problems  have  been  reduced  considerably  from  the  previous  formulations.  These 
reductions  reflect  a  substantial  Improvement  in  our  understanding  of  the  character¬ 
istics  and  capabilities  of  the  optimization  program  and  result  in  a  comparable 
reduction  in  the  effort  required  for  this  process . 

A.  Transmission  Line  Problem 

The  transmission  line  model  has  been  revised  and  simplified  by  the 
elimination  of  all  equality  constraints .  The  parametric  constraint  on  power  has 
been  replaced  by  an  input  constant  -  -  and  the  remaining  constraint  by  an 

auxiliary  equation.  In  addition,  the  voltage  has  been  eliminated  as  a  variable 
and  replaced  by  a  constant  of  100,000  volts.  The  net  effect  has  been  a  reduction 
of  the  model  from  a  system  of  7  Independent  variables  with  4  inequality  constraints 
(each  equality  constraint  was  expressed  as  a  pair  of  inequality  constraints)  to  a 
system  of  4  independent  variables  with  no  constraints.  The  resulting  model  is 
illustrated  in  Table  1 . 
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The  above  revised  model  was  reprocessed  by  the  optimization  program 
in  order  to  correct  the  results  provided  in  the  previous  report  which  had  been 
obtained  with  several  of  the  parameter  interaction  terms  omitted.  The  results 
obtained  are  illustrated  in  Table  2.  These  results  indicate  somewhat  larger 
transmission  line  lengths  and  correspondingly  higher  line  loss  factors,  but  do 
not  differ  appreciably  from  the  results  obtained  with  the  previous  formulation. 

Influence  coefficients  have  been  obtained  to  indicate  the  sensitivity  of 
the  weight  to  changes  in  the  design  parameters  from  their  optimum  values .  The 
Influence  coefficients  have  been  defined  as  the  increase  in  total  system  weight 
resulting  from  a  1%  Increase  in  each  of  the  design  variables.  Results  are  con¬ 
tained  in  Table  3.  The  negative  signs  for  the  temperature,  frequency  and  vol¬ 
tage  Influence  coefficients  indicate  that  further  reductions  in  system  weight 
can  be  obtained  by  increasing  these  parameters  beyond  their  maximum  values . 

If  these  maximum  values  reflect  a  state-of-the-art  limitation,  as  is  the  case 
with  the  temperature  limit  for  example,  the  influence  coefficient  provides  a 
measure  of  the  return  to  be  achieved  from  a  development  program  structured  to 
raise  this. limit.. 

B.  Radiator  Panel  Problem 

The  radiator  panel  model  has  been  similarly  revised  and  simplified.  The 
equality  constraint  equations  have  been  removed  an.d  replaced  by  equivalent 
auxiliary  equations  and  the  variables  defined  by  these  equations  eliminated- 
from  the  list  of  Independent  parameters.  The  net  effect  pf  this  process  is  that 
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bounds  on  these  auxiliary  variables  are  removed.  Although  this  did  not  Introduce 
any  complications  in  the  transmission  line  problem,  the  radiator  optimization 
could  not  achieve  a  realistic  design  without  some  additional  constraints.  It  was 
necessary,  therefore,  to  identify  and  use  these  inequality  constraints  to  replace 
the  constraints  applied  by  the  former  bounds  on  the  auxiliary  parameters. 

The  transform  for  normalizing  temperature  to  the  0  to  1  range  was  retained. 
The  transforms  for  the  remaining  independent  variables  were,  however,  replaced 
by  a  simple  multiplier  which  would  set  the  maximum  value  of  each  normalized  i. 
variable  to  1  but  allow  the  minimum  value  to  float.  The  net  effect  of  this  approach 
is  to  convert  the  standard  transform  equation: 

X  5=  aX  +  b 

to  the  form:  x  =  aX 

which  results  in  a  substantial  reduction  in  the  amount  of  manual  calculations 
required  for  converting  the  results  of  the  model  development  analysis  to  the  form 
required  for  the  optimization  program. 

The  net  effect  has  been  the  reduction  of  the  radiator  model  from  a  system 
of  10  Independent  variables  with  10  inequality  constraints  (5  pairs  of  equality 
constraints)  to  a  system  of  3  input  parameters,  4  independent  variables,  and  3 
inequality  constraints.  The  resulting  model  is  illustrated  in  Table  4. 

Six  computer  optimization  runs  were  made  with  this  model  in  order  to  ex¬ 
plore  the  characteristics  of  the  radiator  panel.  The  results  of  these  runs  are 
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summarized  in  Table  5 .  Runs  1  and  2  were  made  with  a  panel  heat  capacity 
of  150  KW  thermal  and  a  fin  effectiveness  of  70%  but  with  different  Initial 
temperatures .  The  results  would  seem  to  Indicate  that  the  temperature  effect 
Is  not  strong  enough  to  Indicate  a  single  optimum  design.  The  probability  of 
achieving  10,000  hours  of  operation  without  a  single  meteorite  puncture  Is  also 
indicated.  Note  that  this  survival  probability  increases  with  increased  tem¬ 
perature  and  that  the  ambiguity  with  respect  to  temperature  would  have  been 
eliminated  If  the  design  had  been  constrained  to  a  single  value  of  survival 
probability.  This,  however,  was  not  done  In  order  to  simplify  the  Illustrative 
problem . 

A  comparison  of  runs  2,3,  and  4  indicate  the  trend  in  panel  weight  and 
in  survival  probability  as  the  thermal  heat  capacity  of  the  panel  is  Increased 
with  both  fin  effectiveness  and  panel  temperature  held  constant ,  It  is  apparent 
that  the  temperature  should  be  Increased  with  Increased  heat  capacity  in  order 
to  maintain  a  given  level  of  survival  probability. 

Runs  5  and  6  illustrate  the  effects  of  a  variation  in  fin  effectiveness  at 
constant  heat  capacity  and  panel  temperature.  These  results  indicate  a  sub¬ 
stantial  variation  in  panel  weight  but  only  a  relatively  small  effect  on  survival 
probability. 

Although  the  above  data  is  insufficient  to  completely  assess  the  trade¬ 
off's  between  panel  temperature  and  fin  effectiveness  it  is  quite  apparent  that 
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these  two  parameters  are  the  most  significant  variables  in  establishing  an 
optimum  panel  weight  at  a  fixed  level  of  survival  probability.  It  is  equally 
apparent  that  this  trade-off  could  be  completely  resolved  by  either  the  cal¬ 
culation  of  three  to  six  runs  or  by  the  augmentation  of  the  radiator  model  by 
an  additional  constraint  on  the  survival  probability. 

Table  6  contains  a  summary  of  Influence  coefficients  for  each  of  the 
six  computer  runs.  These  were  obtained  by  increasing  the  normalized  values 
of  each  of  the  Independent  variables  by  1%.  The  actual  change  Involved  is 
Indicated  in  parentheses.  Note  that  all  of  the  Influence  coefficients  are 
negative  indicating  that  an  Increase  in  any  of  the  independent  variables  will 
result  in  a  decrease  in  radiator  panel  weight.  The  tube  spacing  would  appear 
to  be  the  most  sensitive  design  variable  over  the  range  of  variation  investigated. 

All  of  the  six  runs  were  forced  up  against  the  third  of  the  three  con¬ 
straints.  This  limit  is  imposed  by  a  requirement  for  an  armour  thickness 
greater  than  .1  Inches.  The  effect  of  temperature  on  this  constraint  formulation 
was,  however,  ignored  in  order  to  simplify  the  expression.  This  resulted  in 
the  observed  armour  thickness  of  .143  Inches  at  the  high  temperature  levels 
of  runs  2,3,  and  4.  The  fact  that  all  optimization  runs  ended  on  this  constraint 
Indicates  that  the  more  rigorous  expression  should  have  been  Included. 
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III„  LEADER  USER’S  MANUAL 


The  LEADER  User's  Manual  will  contain  a  summary  of  the  procedural 
elements  of  the  LEADER  technique  with  particular  emphasis  on  those  manual 
operations  rejquired  for  the  preparation  of  system  and  sub-system  data  for  the 
model  development  and  optimization  computer  programs  and  for  the  evaluation 
of  the  results  obtained.  It  will  consist  of  two  major  sections  -  one  on  model 
development  procedures  and  the  other  on  system  assembly  and  optimization 
procedures.  The  section  on  model  development  procedures  has  been  completed 
and  is  Included  as  a  supplement  to  this  progress  report.  The  section  on  system 
assembly  and  optimization  is  in  progress  and  is  expected  to  be  completed  with¬ 
in  the  next  few  weeks.  A  brief  discussion  of  the  approach  to  be  used  in  this 
section  is  Included  in  the  following.  The  computer  manual  will  be  re-issued 
as  one  of  the  technical  summary  reports  summarizing  the  work  performed  under 
this  study  contract . 

The  first  step  in  the  procedure  is  the  Identification  of  the  objective 
function  -  the  parameter  to  be  optimized  -  and  of  the  over-all  system  constraints 
to  be  imposed.  The  Individual  sub-system  models  will  then  be  analyzed  in 
order  to  identify  those  quantities  required  as  Inputs  to  each  model  and  those 
quantities  that  will  be  available  as  outputs.  The  sequence  of  model  calculations 
must  be  established  in  order  to  ensure  that  sub-systems  are  not  calculated  until 
the  required  inputs  are  available. 
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Relationships  which  link,  two  or  more  subroutines  must  then  be  added. 

This  will  Include  such  considerations  as  heat  balances ,  energy  balances ,  flow 
continuity,  speed  matching,  etc.  Each  of  these  relationships  that  apply  can 
be  used  to  eliminate  one  Independent  variable .  The  variables  to  be  eliminated 
should  be  selected  with  care  since  they  will  become  auxiliary  variablea  which 
cannot  be  maintained  within  specified  boundaries  by  the  optimization  process 
unless  specific  constraint  equations  are  provided. 

All  system  and  sub-system  constraints  must  be  identified.  Equality 
constraints  must  be  treated  in  a  fashion  similar  to  the  equations  defining  the 
auxiliary  variables  and  will,  consequently,  permit  the  elimination  of  an  additional 
Independent  variable .  Inequality  constraints  may  be  imposed  but  should  be 
expressed  as  functions  of  the  remaining  independent  variables.  The  ranges  of 
Interest  of  the  independent  variables  must  be  identified  and  the  variables  normalized 
to  a  common  range .  The  model  and  constraint  equations  must  be  modified  to 
receive  the  normalized  variables. 

The  resulting  equations  must  then  be  prepared  for  Insertion  into  the  opti¬ 
mization  program.  All  expressions  which  deviate  from  the  standard  polynomial 
format  will  involve  special  handling  in  order  to  get  them  into  the  optimization 
program.  Finally,  a  set  of  starting  values  must  be  calculated  which  will  satisfy 
all  of  the  system  constraints  imposed. 
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The  system  can  then  be  cuji  on  the  computer  and  the  optimum  design 
obtained.  The  computer  results  will.  In  addition,  Indicate  the  particular  con¬ 
straints  which  are  In  operation  at  the  resulting  optimum  design  and  the  sensitivity 
of  the  objective  function  to  small  deviations  of  the  independent  variables  from 
the  optimum  design. 

The  above  procedures  are  being  described  In  considerably  more  detail 
In  the  second  section  of  the  User's  Manual  which  is  In  preparation.  These 
descriptions  are,  in  addition,  being  Implemented  by  examples  which  will  aid  In 
Illustrating  each  specific  operation.  This  section  will  then  be  combined  with 
the  existing  section  on  model  development  procedures  and  the  report  re-lssued 
as  a  technical  summary  report.  The  computer  program  listings  will  be  Included 
as  appendices  to  this  report  In  order  to  facilitate  the  utilization  of  the  over¬ 
all  technique . 
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IV.  SPUR  MODEL  FORMULATION 


Effort  has  been  continued  on  the  formulation  of  the  350  KW  arc  jet  and  the 
1  MW  ion  engine  power-propulsion  system  of  SPUR  during  this  reporting  period. 
These  efforts  have  been  associated  previously  with  the  accumulation  of  basic 
design  data  and  performance  characteristics,  with  the  development  of  sub-system 
design  procedures,  and  with  the  identification  of  appropriate  design  criteria  and 
assumptions.  This  process  has  been  completed  for  the  turbine  and  the  generator 
sub-systems  and  is  expected  to  be  completed  for  the  remaining  sub-systems 
within  the  next  few  weeks. 

The  results  of  these  analyses  have  been  used  to  develop  preliminary  model 
formulations  for  the  turbine  and  generator  sub- systems.  These  formulations  serve 
as  a  summary  of  the  results  of  the  respective  design  analyses  and  as  an  indication 
of  the  additional  analysis  required  to  reduce  the  design  procedure  to  the  form 
required  by  the  optimization  program.  The  preliminary  turbine  and  generator 
formulations  will  be  described  in  the  following  section  along  with  the  design 
approaches  and  assumptions  upon  which  they  are  based. 

The  overall  approach  has  been  based  upon  the  identification  of  sub-system 
weight  as  the  objective  function.  The  Individual  sub- system  weights  will  then 
be  summed  and  the  overall  weight  minimized  by  the  optimization  process.  This 
minimization  will  be  carried  out  subject  to  constraints  on  overall  system  length 
and  diameter  required  for  integration  with  the  Saturn  C-IB  boost  vehicle. 
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A.  TURBI^E  SUB-SYSTEM 


The  turbine  design  philosophy  had  been  based  upon  the  use  of  conventional 
steam  turbine  practice  employing  full  admission-root  impulse  stages.  The  follow¬ 
ing  design  assumptions  have  been  made: 

1.  The  flow  of  saturated  vapor  can  be  represented  with  sufficient  accuracy 
by  the  use  of  perfect  gas  formulae  when  suitable  values  of  the  gas 
constant  and  the  specific  heat  ratio  are  employed. 

2.  The  total  to  total  efficiency  of  a  wet  vapor  turbine  can  be  represented 
as  a  function  of  the  average  stage  pitch  line  velocity  ratio  Up/Cp 
(where  Up  in  the  pitch  line  velocity  and  Cq  is  the  velocity  associated 
with  converting  the  total  enthalpy  drop  across  the  stage  to  a  kinetic 
energy),  the  interstage  moisture  removal  effectiveness,  and  the  tip 
clearance  to  bucket  span  ratio.  The  data  of  Figure  1  illustrates  the 
effect  of  the  velocity  ratio  and  the  moisture  removal  factor.  The 
efficiency  obtained  from  Figure  1  is  corrected  by  subtracting  the 
quantity: 

3  (Clearance  to  Span  Ratio) 

This  correlation  assumes  a  high  level  of  design  and  manufacturing 
precision  on  bucket  profiles  and  interstage  labyrinth  seals  and  has 
been  obtained  from  a  large  body  of  steam  turbine  test  data. 

3.  Turbine  discs  are  of  the  constant  stress  type. 
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4. '  Refractory  metal  alloys  are  used  in  turbine  rotor  construction  which 

permit  the  use  of  design  stress  limits  for  0,1%  creep  after  10,000 
hour  operation  as  shown  by  the  dSta  of  Figure  2. 

5.  Untapered  buckets  which  are  integral  with  the  wheel  have  been  used. 

A  schematic  drawing  of  the  turbine  is  included  in  Figure  3.  Table  7  contains 
a  summary  of  the  preliminary  turbine  model  developed  which  has  been  arranged 
according  to  the  following  format: 

1.  Input  -  The  independent  variables- which  establish  the  turbine  design 
in  sufficient  detail  to  permit  an  estimation  of  its  weight.  These 
parameters  will,  in  general,  be  optimized  by  the  LEADER  technique 
unless  they  are  obtained  as  auxiliary  outputs  from  the  other  sub- system 
models. 

2.  Auxiliary  Equations  -  The  dependent  variables  which  are  needed  for  the 
weight  calculation,  for  establishing  the  sub- system  constraints,  or  as 
inputs  to  the  other  sub-system  models. 

3.  Constraint  Equations  -  The  physidal  constraints  imposed  by  turbine 
design  considerations.  These  constraints  must  be  augmented  by  linear 
constraints  and  by  additional  non-linear  constraints.  The  linear  con¬ 
straints  are  the  bounds  that  are  established  on  the  Independent  variables. 
The  non-linear  constraints  are  introduced,  where  needed,  in  order  to 
establish  bounds  on  the  dependent  variables.  These  must  be  expressed 
in  terms  of  two  or  more  independent  variables. 
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4.  Objective  Function  -  A  summary  of  the  Individual  turbine  element 
weights  which  will  be  added  to  the  weights  of  the  other  sub-systems 
and  minimized. 

5.  Parameters  -  Those  factors  which' will  be  assumed  to  be  constant 
during  any  optimization  run.  Any  of  these  parameters  could  be  varied 
in  subesequent  optimization  runs  or  transferred  to  the  input  section 
where  they  would  be  treated  as  additional  independent  variables. 

6.  Empirical  Relations  -  Empirical  equations  required  for  the  evaluation 
of  the  auxiliary  variables  which  must  be  developed  by  the  techniques 
described  in  the  LEADER  User's  Manual.  Note  that  the  turbine  equa¬ 
tions  required  are  functions  of  either  one  or  two  variables  for  each 
working  fluid  selected. 

7.  Auxiliary  Outputs  -  Parameters  which  are  expected  to  be  required  as 
input  to  other  sub-system  models.  This  list  will  be  augmented  as  the 
other  model  formulations  are  developed. 

The  above  format  has  been  utilized  as  ah  aid  in  Identifying  the  interfaces 
between  the  turbine  sub-system  and  the  other  sub-systems  and  in  identifying 
the  specific  variables  which  will  be  optimized  by  the  LEADER  technique.  It  is 
particularly  important  to  keep  a  running  count  on  the  number  of  independent  variables 
as  a  total  in  excess  of  50  would  require  the  overall  system  model  to  be  treated  as 
two  or  more  sub-optimization  problems. 
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B.  GENERATOR  SUB-SYSTEM 


The  generator  sub-system  model  has  been  developed  for  a  radial  gap  generator. 
A  similar  model  is  in  process  for  an  axial  gap  generator.  The  radial  gap  design 
procedure  has  been  based  upon  the  following  assumptions: 

1.  The  electromagnetic  design  is  limited  by  values  of: 

15  for  the  ratio  of  rotor  tooth  depth  to  flux  gap  length  and 
30  for  the  ratio  of  rotor  tooth  gap  at  the  outside  of  the  tooth 
to  flux  gap  length. 

2.  Corresponding  to  the  above  limits,  the  stator  to  rotor  leakage  flux 
has  been  assumed  to  be  20%  of  the  working  gap  flux. 

3.  The  effective  A.  C.  flux  is  1/2.2  times  the  working  gap  flux. 

4.  The  pole  embrace  is  2/3. 

5.  Rotor  hot  spots  are  within  50°  of  the  coolant  temperature,  and,  conse¬ 
quently,  the  rotor  stress  limit  can  be  expressed  as  a  function  of  the 
coolant  outlet  temperature.  Assumed  values  are  illustrated  in  Figure  4. 

6.  The  field  coil  proportions  have  been  selected  to  correspond  with  mini¬ 
mum  frame  weight  per  coil  cross  sectional  area  and  mean  diameter. 

7.  Stator  stack  and  conductor  weights  have  been  calculated  by  applying 
a  correction  factor  to  account  for  the  end  turns  and  end  turn  potting 
to  an  armature  weight  based  upon  a  solid  lamination  material 

A  schematic  drawing  of  the  generator  is  included  in  Figure  5.  Table  8  contains 
a  summary  of  the  preliminary  generator  model  developed  according  to  the  same 
format  used  for  the  turbine. 
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TABLE  1 


TABLE  1  (Continued) 


Transmission  Line  Modal 


Objective  Function 


W  =  12,277.2  +  19,135.5  Pi  +  4741.5  Pixi  -2423.6  PiX2  -  1846.3  Pjxs 

+  689.2  P1X4  -599.9  PjXiX2  -  457.0  P1X1X3 
+  170.6  P1X1X4  +13,094.4  P1X4G1 
+  771.9  xi  -  97.7  X1X2  -  74.4  xjxg  +  27.8  XIX4  +  527.6  X1X4G1 

-  769.5  X2 

-  3400 . 6  X3 

+  112.2x4  +  2131.6  X4G1 

-  7635.4  G2  +  2306.4  G2G3  -  69.2  G2G4  +  2616.6  G2^ 

-  1495.8  G25G3  +  448.7  G2^G4 

+  4317.4  G3 
+  129.5  G4 
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TABLE  2 

•j-f 

!! 

OPTIMUM  OVER-ALL  SYSTEM  DESIGN 


Transmission  Line  Model 


Power  -  P 

KW 

350 

1000 

2500 

Efficiency  -  n 

% 

94.20 

96.48 

93.74 

Frequency  -  f 

cps 

3199.8 

3055.2 

3058.6 

Temperature  -  T 

2500.0 

2500.0 

2500.0 

Length  -  1 

ft 

71,880 

44,090 

50,100 

Voltage  -  E 

Volts 

lOlO 

o 

o 

lOlO 

Loss  Factor  -  l-ri 

% 

5.80 

3.52 

6.26 

Total  System  Weight  - 

W  lbs 

4,845.3 

10,179 

22,174 

Specific  Weight  -  W/P 

Ibs/KW 

13.84 

10.18 

8.86 

TABLE  3 


INFLUENCE  COEFFICIENTS 


Transmission  Line  Model 


Power  Output 

KW 

350 

1000 

Total  System  Weight 

lbs 

4,845.3 

10,179 

Influence  Coefficients 

-  lbs 

Power 

168.6 

169.1 

Efficiency 

.111 

67.50 

Temperature 

-36.36 

-39.89 

Frequency 

-9.41 

-19.42 

Line  Length 

-1.90 

.8000 

2500 

22,175 


170.3 

123.5 

>53.88 

-32.09 

1.00 


Voltage 


-1.572 


-1.733 


6.34022 
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TABLE  5 


Optimum  Radiator  Panel  Design 


Run  No. 

1 

2 

3 

4 

5 

6 

Heat  Capacity 

0 

KWT 

150 

150 

300 

450 

450 

450 

Fin  Effectiveness 

Tli 

.70 

.70 

.70 

.70 

.80 

.90 

Heat  Transfer  Coefficient 

2 

.41 

.41 

.41 

.41 

.20 

.08 

Temperature 

T 

°R 

1671 

1754 

^754 

1753 

1678 

1678 

Fin  Thickness 

t 

in 

.120 

.148 

.149 

.150 

.114 

.114 

Tube  Diameter 

d 

in 

.376 

.256 

.253 

.252 

.382 

.382 

Tube  Spacing 

N/L  ft 

8.35 

8.49 

8.49 

8.48 

10.28 

12.59 

Vapor  Volume  Flow 

V 

cfs 

14.76 

10.26 

20.62 

31.02 

43.24 

43.17 

Number  of  Tubes 

N 

31.9 

48.2 

98.6 

149.9 

90.9 

90.4 

Fin  Length 

1 

In 

.422 

.435 

.438 

.439 

.285 

.180 

Armour  Thickness 

6 

In 

.108 

.144 

.143 

.143 

.107 

.105 

Header  Diameter 

D 

In 

3.68 

3.07 

4.54 

5.34 

6.30 

6.29 

Panel  Length 

L 

ft 

3.82 

5.65 

11,62 

17.65 

8.84 

7.18 

Panel  Width 

h 

ft 

6.36 

3.46 

3.24 

3.10 

7.29 

8.53 

Panel  Area 

A 

ft2 

48.6 

39.0 

75.2 

109.4 

128.8 

122.4 

Panel  Vulnerable  Area 

Av 

ft2 

23.6 

15.7 

35.0 

55.3 

80.8 

88.9 

Survival  Probability  In 

P 

% 

87.3 

96.1 

91.4 

86.8 

62.8 

59.9 

10,000  hrs. 

W  lbs  42.4  43.2  94.9  151.8  127.1  118.6 


Radiator  Panel  Weight 


TABLE  6 


INFLUENCE  COEFFICIENTS 
Radiator  Panel  Model 


Run  No. 

1 

2 

3 

m 

5 

6 

Panel  Heat  Capacity 

-  KWT 

150 

150 

300 

450 

450 

450 

Fin  Effectiveness 

.70 

.70 

.70 

.70 

.80 

.90 

Panel  Temperature  - 

«R 

1671 

1754 

1754 

1753 

1687 

1687 

Influence  Coefficients  -  lbs 

Fin  Thickness  -  In 

(A  =  ,003  In) 

-.0022 

-.0124 

-.015 

-.0035 

-.066 

-.048 

Tube  Spacing  -  ft"^ 

(A  *  .67  ft’^) 

-.56 

-.64 

-1.28 

-2.08 

-1.70 

-1.35 

Tube  Diameter  -  In 

(A  =  .005  In) 

-.057 

-.040 

-.171 

-.024 

-.140 

-.137 

Temperature  -  ®  R 

(A  =  2.1°R) 

-.024 

1 

o 

o 

00 

-.081 

-.0096 

-.079 

-.060 

TABLE  7 


Input 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 

Auxiliary 

1. 

2. 

3. 

4. 

5. 

6. 

7. 


Preliminary  Turbine  Sub-Svstem  Model 


Turbine  Inlet  Total  Temperature  -  T^^  -  °  R 
Turbine  Inlet  Total  Pressure  “ 

Mass  Flow  Rate  -  W  -  Ib/sec 
Nozzle  Angle  -  a 

First  Stage  Bucket  Tip  Radius  -  r^^  -  in 
First  Stage  Nozzle  Exit  Velocity  ~  -  fps 

Leaving  Loss  Factor  -  LL 
Moisture  Removal  Parameter  -  M 
Clearance  -  CL  -  in 

Turbine  Discharge  Total  Temperature  "  ^ 


Equations 


‘hi 


N 


bl 


»/7g1R5, 


tl 


FF^  =  f 


i^i] 


*1  = 


Pti  (PFj)  sin  a 

"w/kT 


-yT 


"1 


30  U 


hi 


"'^hl^tl 


4.52 


/ 


N 

irnju, 


(First  Stage  Nozzle  Velocity  Function) 
(First  Stage  Flow  Function) 

(First  Stage  Annulus  Area) 

(First  Stage  Rotor  Hub  Radius) 

(First  Stage  Hub  Velocity) 

(Rotor  Speed) 

(First  Stage  Bucket  Stress) 


TABLE  7  (Continued) 


Preliminary  Turbine  Sub-System  Model 
B.  Auxiliary  Equations  (Continued) 


8. 

F 

m 

— 

[Working  Fluid,  T^^,  T^^  3 

(Ratio  of  Ideal  Ranklne  to 
Carnot  Efficiency) 

9. 

H 

V 

s 

[Working  Fluid,  T^^  ] 

(Heat  of  Vaporization) 

10. 

^t 

s 

.85  H  F  [1-T_/T,  ] 

(Nominal  Enthalpy  Drop) 

11. 

92 

= 

/  2gJA.H^LL 

(Leaving  Velocity) 

12. 

^2 

= 

Cz/  y 

(Last  Stage  Velocity  Function) 

13. 

^2 

= 

^4  ^  ^t2  ^ 

(Total  Discharge  Pressure) 

14. 

FF2 

s 

fi  [v^  ] 

(Last  Stage  Flow  Function) 

15. 

^2 

= 

[P^2  (FF2)]/[W^T-~] 

(Last  Stage  Annulus  Area) 

16. 

^h2 

= 

.5C2  /  sin  a 

(Last  Stage  Hub  Velocity) 

17. 

'^h2 

= 

[30  1/  [/"AjN^  *900  U 

2 

j^2  ]  (Last  Stage  Hub  Radius) 

18. 

^2 

= 

’'>"h2 

"^'•h2 

(Last  Stage  Tip  Radius) 

19. 

.5  gjAH, 

20. 

n 

= 

[2IU^']/[Uhi'.\2'] 

(No.  of  Stages) 

21. 

L 

2 .5  n  S 

(Rotor  Length) 

22. 

"l 

= 

^5  ^"^tl^ 

(First  Stage  Allowable  Stress) 
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TABLE  7  (Continued) 
Preliminary  Turbine  Sub-System  Model 

B .  Auxiliary  Equations  (Continued) 


23. 

‘^nl 

(First  Stage  Wheel  Neck 
Thickness) 

24. 

‘hi 

= 

nl 

(First  Stage  Wheel  Hub 
Thickness) 

25. 

= 

4.52A  |.  (N/IOOO)^ 

(Last  Stage  Bucket  Stress) 

26. 

's' ^<2* 

(Last  Stage  Allowable  Stress) 

27. 

‘n2 

= 

[=rb2S]/[2F2] 

(Last  Stage  Wheel  Neck 
Thickness) 

28. 

‘h2 

= 

t,e^»b«hVr^  =  -b2'’^ 

n2 

(Last  Stage  Wheel  Hub 
Thickness) 

29. 

= 

(First  Stage  Pitch  Line 
Velocity  Ratio) 

30. 

CLR 

= 

o 

f 

1 _ 1 

\ 

t — 1 

t 

1 _ 1 

(First  Stage  Clearance  Ratio) 

31. 

nti 

= 

fg  [Vj,  M  ]  -  3  (CLR)^ 

(First  Stage  Efficiency) 

32. 

= 

(Last  Stage  Pitch  Line 
Velocity  Ratio) 

33. 

CLR^ 

= 

rCL]/[rt2^rh2] 

(Last  Stage  Clearance  Ratio) 

34. 

^t2 

= 

fg  [V2,  M]-  3(CLR)^ 

(Last  Stage  Efficlenay) 

35. 

^t 

s 

[’'i''hr''2<2]/t4 

(Over-all  Turbine  Efficiency) 
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TABLE  7  tSontlnued) 


Preliminary  Turbine  Sub-System  Model 


B.  Auxiliary  Equations  (Continued) 


36. 

P 

= 

i:(Ti^/.8$  -  LL] 

(Turbine  Power  Output) 

37. 

^bl 

= 

•5  (4  -  4) 

(First  Stage  Bucket  Weight) 

38. 

^1 

z= 

2/3 

(First  Stage  Rim  Weight) 

39. 

W  , 
wl 

= 

(First  Stage  Wheel  Weight) 

40. 

W  , 
nl 

= 

(First  Stage  Nozzle  Weight) 

41. 

=; 

[4npstg/TrA^ 

(Inlet  Satoll  Weight) 

42. 

a 

w. ,  ^  w  ,  w  ,  W  , 

bl  rl  wl  nl 

(First  Stage  Weight) 

43. 

W 

b2 

= 

•  5  Hf^S  (r^2  ■  ^h2^ 

(Last  Stage  Bucket  Weight) 

44. 

W 

^r2 

= 

2/3 

(Last  Stage  Rim  Weight) 

45. 

W  „ 
w2 

= 

'b'^h2  ^‘h2  *  ^  *n2 

(Last  Stage  Wheel  Weight) 

46. 

W 

^32 

a 

(Exit  Scroll  Weight) 

47, 

^n2 

= 

2»s'*'b2/^ 

(Last  Stage  Nozzle  Weight) 

48. 

W2 

= 

'^b2  *^2  '«'w2  *'^n2 

(Last  Stage  Weight) 

49. 

W 

c 

= 

3/4  nP3t^(rjj  ‘  *  tj  L 

(Casing  Weight) 

Constraints 

1. 

a  <F 
bl  1 
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TABLE  7  (Continued) 


Preliminary  Turbine  Sub-Svstem  Model 

C .  Constraints  (Continued) 

D.  Objective 

1.  W  =  2  Wgj  ^  n/2  (Wj  ^  W^)  ^  2  Wg2  ^ 

E.  Parameters 

Gas  Constant  -  R 

Density  of  Rotor  Material  -  (Molybdenum) 
Density  of  Stator  Material  -  p  (Columbium) 
Bucket  Chord  -  S 
Scroll  Thickness  -  t 

s 

Casing  Thickness  -  t^ 

Specific  Heat  Ratio  -  y,  1.22  for  Potassium 

F.  Emperical  Relations 


1. 

FF 

fj  (V) 

2. 

F 

m 

(Working  Fluid,  T^^,  T^2) 

3. 

H 

V 

f^  (Working  Fluid,  T^^) 

4. 

^4  ^^t2^ 

5. 

F 

's'V 

6. 

Tit  " 

fg  (Vj,  M) 

G.  Auxiliary  Outputs 

Rotational  Speed  -  N,  RPM 
Turbine  Power  Output  -  P,  BTU/  sec 
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TABLE  8 


Preliminary  Generator  Sub-Systen 
(Radial  Gap  Machine) 


A.  Input 


1. 

Rotational  Speed  -  N ,  RPM 

2. 

Frequency  -  £,  cps 

3. 

Voltage  -  E 

4. 

Gap  Diameter  -  D  ,  in 

g 

5. 

Gap  Length  -  L  ,  in 
g 

6. 

Flux  Density  Across  Gap  - 

B 

g 

7. 

Depth  of  Core  Behind  Slots 

-  d,  in 

8. 

Slot  Depth  -  S,  in 

9. 

Slot  Pitch  -  p,  in 

10. 

Current  Density  -  ID 

11. 

Load  Power  Factor  -  PF 

12. 

Coolant  Temperature  -  T  , 

c 

°R 

13. 

Deep  Bar  Factor  -  F^ 

14. 

Turbine  Power  Output  -  P , 

watts 

Auxiliary  Eauatlons 

1. 

B,-  =  .4nD  L  B 

tf  g  g  g 

2. 

] 

Model 


(Frame  Flux) 
(Core  Diameter) 
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TABLE  8  (Continued) 


Preliminary  Generator  Sub-System  Model 
(Radial  Gap  Machine) 

B.  Auxiliary  Equations  (Continued) 


3. 

g 

= 

[D-D  ]/30 
‘  g  c 

(Gap) 

4. 

T 

= 

60  f/N 

(No.  of  Rotor  Teeth) 

5. 

s 

s 

nD^[6  Tp] 

(Slots/Pole/Phase) 

6. 

I 

= 

.9  P/[E(PF)  ] 

(Current) 

7. 

n 

= 

[4  (10)‘®fSK^KpB^^]/  [1.2  E  ] 

(No.  of  Parallel  Paths 
thru  Winding) 

8. 

IS 

= 

2I/n 

(Current  per  Slot) 

9. 

w 

s 

[IS]  /  [  ID(SSF)  S] 

(Slot  Width) 

10. 

"st 

z: 

[BgP  (FF)]/  [p  -  w] 

(Flux  Density  in  Stator  Teeth) 

11. 

®sc 

= 

B^^  /[2.4Td] 

(Flux  Density  in  Stator  Core) 

12. 

AT 

= 

.626  F^gBg  ^  1.35  (IS) 

(Ampere  Turns) 

13. 

‘f 

= 

[2  (AT)  (  FSF)  /  (IDD) 

(Field  Coil  Length) 

14. 

r 

w 

tt/4  PgWlOOO./ [1.507  (Dg  - 

)  D  /D^  *  1.230  ttD^] 
g  c  c-' 

(Rotor  Stress) 

15. 

P 

= 

[T  ] 

1 c 

(Resistivity) 

16. 

R 

s 

s 

[  pF^(2  Lg  +  .55  n  D^/T]  /  Wg(S  SF ) 

(Resistance  per  Slot) 
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TABLE  8  (Continued) 


Preliminary  Generator  Sub- System  Model 
(Radial  Gap  Machine) 


17.  =  6TSRg(IS)^  (Conductor  I^R  loss) 

18,  P,  =  2np^L  K„(FLT)2{sDgB3,2[_E_|.  ,2583^^ 

P“W  " 

[(Dg  +  +  2d)^  -  (Dg  +  (Stator/lon  Loss) 


19. 

KPc 

>  Pg ,  go  to  line  21 

20. 

V 

=  Pg,  go  to  line  22 

21. 

Ppf 

=  Pc 

(Rotor  Pole  Face  &  Stray  Load  Loss) 

22. 

R,= 

2np^[Dg  t  2s  t  2d  +  .  51f]/[(FSF)lf 

2 

]  (Field  Resistance) 

23. 

Pf 

-  (AD^Rf 

(Field  Excitation  Loss) 

24. 

PG 

=  P-  [Pg  +  Pg  Ppf  i.  Pf] 

(Generator  Power  Output) 

25. 

EFF 

=  PG  /  P 

(Generator  Efficiency) 

26. 

^fr 

=  (21f  ^  Lg)p3  Bff  /  Bf 

(Frame  Weight) 

27. 

^fc 

TJ  /  2  Pc^f^  2g  -f  2d 

.51f] 

(Field  Coll  Weight) 

28.  Wg  =  .SSTTPgLg  [(Dg  4  2s  +  2d)2  -  Dg2] 

(Staok  Weight) 

29.  Wr  =  rr  /  4  pg  [Dg2  (2lg  4  If)  ^  2/3  Lg(Dg2  - 

(Rotor  Weight) 
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TABLE  8  Continued 


C.  Constraints 


Preliminary  Generator  Sub-System  Model 
(Radial  Gap  Machine) 


1. 

PG  >  350,000 

2. 

2  TT  Dg>  90gT 

3. 

Bst  <Bs) 

4. 

C 

600  Kp  Dg  Bf 

d  (Dg  2d  -f  2s) 


2  ^ 


D.  Objective 

1 .  W  =  Wfr  +  Wf  ^  f  W  +  W 


E.  Parameters 

Density  of  Copper  -  p 

c 

Density  of  Steel  -  p 

s 

Maximum  Flux  Density  across  Rotor  -  (Bj.)^^^ 
Maximum  Flux  Density  across  Stator  -  (85)^1 
D.  C.  Current  Density  across  Field  Coil  -  IDD 
(Frequency)  (Lamination  Thickness)  -  FLT 
Flux  Density  in  Frame  -  Bj 
Fringing  Factor  -  FF 
Slot  Space  Factor  -  SSF 
Field  Space  Factor  -  FSF 
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TABLE  8  (Continued) 


Preliminary  Generator  Sub-Svstem  Model 
Winding  Constant  -  K^Kp  (Radial  Gap  Machine) 

Core  Loss  Proportion  Constant  - 

Ampere  Turn  Factor  -  Fj 

F.  Empirical  Relations 

1.  Resistivity  -  p  =  T^ 

G.  Auxiliary  Outputs 
Current  -  I,  amp. 

Generator  Power  Output  -  PG,  watts 


Generator  Efficiency  -  EFF 


o 
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Figure  2  TURBINE  ROTOR  DESIGN  STRfeSS  LIMITS 
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PRELIMINARY  TURBINE  SUB-SYSTEM  MODEL 
(Multi-Stage  Machine) 


CL  -  Clearance,  in 
Fiji  -  Tip  Radius,  in 
rjj  -  Hub  Radius,  in 
S  -  Chord,  in 


Figure  3  TURBINE  SCHEMA 
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LEADER  USERS  MANUAL 


I.  INTRODUCTION 

The  LEADER  technique  is  the  name  applied  by  the  General  Electric  Company 
to  a  procedure  for  determining  the  optimum  values  of  the  design  variables  of 
highly  complex  engineering  designs.  The  procedure  involves  the  development 
of  an  analytical  model  of  the  complete  system  which  defines  the  interactions 
between  each  of  the  design  variables  and  sub-systems  and  contains  all  desired 
design  constraints  and  physical  system  limitations.  The  resulting  model  Is 
then  explored  numerically,  using  the  method  of  steepest  descent,  and  an 
optimum  design  identified  which  Is  compatible  with  all  of  the  imposed  constraints. 

The  overall  technique  consists  of  two  digital  computer  procedures  -  for 
model  development  and  model  optimization  -  and  supplementary  manual  opera¬ 
tions  for  the  preparation  of  the  Input  data  required  by  the  computer  programs 
and  for  the  evaluation  of  the  results  obtained.  This  manual  contains  a  descrip¬ 
tion  of  these  manual  operations  as  applied  to  a  typical  problem  In  order  to 
0' 

facilitate  the  use  of  the  LEADER  technique.  This  manual  Is  Intended  as  a 
supplement  to  References  (1)  and  (2)  which  describe  the  development  of  each  of 
the  computer  programs,  their  format,  and  their  listings.  Such  detail  has, 
accordingly,  been  omitted  from  this  manual. 

A  glossary  of  terms,  nomenclature,  and  examples  are  Included  in  the 
Appendices  to  this  manual. 
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II.  RESULTS  ACHIEVABLE 


The  results  of  the  model  development  phase  of  the  LEADER  technique 
are  the  reduction  of  the  design  characteristics  of  the  system  under  Investi¬ 
gation  Into  a  completely  analytic  form  through  the  process  of  empirical 
approximation.  The  validity  of  the  approximation  Is  Indicated  through  several 
"goodness  of  fit"  criteria  which  can  be  used  to  substantiate  the  final  results 
of  the  overall  analysis.  In  addition,  the  relative  Importance  of  each  of  the 
variables  and  groupings  of  variables  Is  Indicated  In  order  to  add  some  addi¬ 
tional  Insight  Into  the  characteristics  of  the  overall  system. 

Those  elements  of  the  design  which  can  be  described  analytically,  need 
not  be  represented  unless  their  analytical  expressions  are  too  complex  to  be 
handled  straightforwardly.  Thus,  this  phase  would  be  concerned  primarily 
with  experimental  data,  graphical  data,  or  statistical  data. 

The  results  of  the  system  optimization  phase  are  the  values  of  each  of 
the  Independent  and  dependent  variables  at  the  point  which  optimizes  the  design 
subject  to  the  constraints  specified.  If  the  resulting  design  Is  constrained  by 
any  boundaries  Imposed  upon  the  system,  they  will  be  Identified.  The  sensi¬ 
tivity  of  the  resulting  design  to  small  changes  In  each  of  the  Independent 
design  variables  will  be  Indicated.  Additional  sensitivity  Information  can 
sometimes  be  obtained  from  the  path  followed  by  the  design  In  the  Iterative 
optimization  procedure. 
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III.  MODEL  DEVELOPMENT  PROCEDURES 


The  initial  step  of  the  model  development  process  is  the  collection  of 
all  design  procedures,  experimental  data,  graphical  data,  etc.  which  define 
the  design  problem,  the  interactions  between  all  of  the  significant  parameters, 
and  the  physical  constraints  or  boundaries  imposed  upon  the  system.  In  order 
to  facilitate  the  subsequent  analysis ,  the  system  should  be  broken  down  into 
sub-systems  or  smaller  component  parts.  All  elements  which  can  be  described 
in  analytical  form  should  be  put  aside  and  reserved  for  the  optimization  process. 
All  remaining  non-analytical  elements  must  be  reduced  to  some  type  of  analyti¬ 
cal  form. 

The  procedure  to  be  used  in  developing  the  analytical  forms  for  the 
remaining  elements  will  be  dependent  upon  whether  1,  2,  or  more  independent 
variables  are  Involved.  It  is  desirable,  therefore,  to  subdivide  the  elements 
into  the  smallest  elements  which  will  describe  each  physical  effect  so  that  the 
task  of  empirical  approximation  will  be  minimized.  Each  of  these  procedures 
is  described  in  the  following  sections. 

A.  One  Independent  Variable 

Those  elements  which  can  be  described  in  terms  of  a  single-line 
graph  or  a  single-entry  table  can  be  empirically  represented  by  the  procedure 
contained  in. this  section.  The  result  will  be  an  empirical  equation  of  the  form: 

y  =  Bo  +  Bi  X  f  82  X  ^  f  B3  x^  +  B4  x'^  (1) 
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where  the  polynomial  form  Is  utilized  in  order  to  achieve  results  which  are 


most  easily  compatible  with  the  optimization  procedures  of  the  following 
section. 

The  Initial  step  is  to  select  the  range  of  variation  in  x  over 
which  the  dependent  variable  -  y  -  is  to  be  fitted  and  to  determine  the  degree 
of  approximation  which  will  be  required  in  order  to  achieve  a  satisfactory 
representation.  The  independent  variable  is  normalized  over  the  range  of 
0  to*l  by  the  transformation  equations: 

®  =  Wmln 
=  Wmax 

c  =  — i— 
b-a 

X  =  c(x-a) 

If  the  data  is  in  tabular  form,  it  will  be  necessary  to  plot  it 
before  proceeding.  A  degree  of  polynomial  representation  is  then  selected 
according  to  the  shape  of  the  data.  If  the  data  is  non-linear,  for  example, 
at  least  the  second  degree  should  be  selected.  From  the  data  of  Table  la, 
specific  values  of  X  are  noted  at  which  corresponding  values  of  y  must  be 
obtained  from  its  curve.  Table  lb  is  then  utilized  to  calculate  the  probable  ■ 
fitting  error  associated  with  the  degree  of  polynomial  representation  selected. 

^  ■  ^0^0  ^  ^1^1  ^  ^2^2  ^  ^3^3  ^  •  (3) 


„4_ 


Errors  of  the  order  of  1  to  2  %  are  normally  acceptable.  If  the  calculated 
error  is  considered  unacceptable,  the  process  should  be  repeated  with  the 
next  higher  degree  of  representation.  If  the  use  of  the  fourth  degree  is  still 
unable  to  produce  acceptable  representation  it  will  be  necessary  to  either 
transform  the  variables  or  to  reduce  the  fitting  region. 

Transformations  such  as  Inx,  Iny,  l/x  ,  l/y  ,  sin  x,  sin  y,  etc. 
may  be  utilized  in  place  of  the  original  x,y  values  in  equations  (1),  (2),  and 
(3).  Unfortunately,  no  specific  guidelines  can  be  offered  as  to  the  trans¬ 
formations  that  will  work  for  any  one  problem.  The  user  must  apply  his  own 
judgment  in  this  area.  If  the  transformation  approach  does  not  produce 
acceptable  results,  however,  the  fitting  accuracy  may  be  Improved  by  reducing 
the  range  -  l/c  -  over  which  the  fit  Is  to  be  obtained. 

‘  After  an  acceptable  degree  of  fit  has  been  obtained,  the  data  of 

Table  2  Is  used  to  obtain  the  actual  coefficients  of  the  fit.  For  a  second 
degree  fit ,  for  example ,  the  first  coefficient  A©  is  obtained  from  the  equation : 

Aq  =  ,  8333 . .  •  y^  +  .  333.  ••  y^  ~  .  333. f  .  166. .  .y^  (4) 

Similarly,  the  values  of  A^  and  A2  are  obtained  from  the  other  lines  of  the 
'second  degree  table.  The  resulting  equation  Will  express  the  dependent 
variable  in  terms  of  the  normalized  Independent  variable: 

y  =  Ao+AiX  +  A2X^  +  . ..  AnX"  (5) 
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Equation  (5)  must  then  be  unnormalized  to  reduce  it  to  the  form 
of  equation  (1).  This  can  be  accomplished  by  calculating  columns  2  thrdugh 
9  of  Table  3.  The  first  unnormalized  coefficient  is  then  obtained  from  a  sum¬ 
mation  of  the  products  of  column  4  with  column  5: 

=  Aq  -  acAj  +  -  a^c^a3  +  a^c^A^  (6) 

Similarly,  the  other  coefficients  are  obtained  from  the  product  of  column  4 
with  each  of  the  remaining  columns  6  through  9. 

The  resulting  equation  can  then  be  used  directly  in  the  opti¬ 
mization  process.  The  curve  fitting  technique  employed  minimizes  the  amount 
of  manual  effort  required  for  obtaining  the  required  coefficients.  It  is 
described  in  more  detail  in  Reference  (3).  If,  however,  computer  curve  fitting 
programs  are  available ,  they  may  be  used  in  place  of  the  above  procedure . 

B.  Two  Independent  Variables 

Those  elements  of  the  system  w4iich  can  be  described  in  terms 
of  a  family  of  curves  or  a  double-entry  table  can  be  represented  by  a  proce¬ 
dure  which  is  an  extension  of  the  one  described  in  Section  A.  The  result  will 
be  an  empirical  equation  of  the  form: 

y  =  Boo  +  ^01  ®02  +  •••  Bon^'' 

+  BiqZ  +  Bii  xz  +  Bi2X^z  +  . . .  'in*"'  (7) 

+  *20'^  *  '21**^  ♦  ®22*^*  ♦  •  •  •  '2n**'^ 

+  -  «. 
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The  initial  step  is  to  determine  the  range  of  variation  required 
in  X  and  z  and  to  determine  the  degree  of  approximation  required  for  each 
variable.  Note  that  the  degree  in  x  need  not  be  the  same  as  the  degree  in  z. 
Both  variables  must  be  normalized  over  the  0  to  1  range  by  equation  (2). 

Table  1  is  then  used  to  determine  fitting  errors  In  both  the 
X  and  z  directions.  Transformation  techniques  must  be  considered  at  this 
point  if  desirable  fitting  accuracies  are  not  obtained. 

After  the  required  degree  of  fit  in  each  direction  has  been 
determined ,  the  coefficients  are  determined  in  a  manner  quite  similar  to 
that  used  in  the  unl-varlate  fitting  process  described  in  the  previous  section. 
The  bi-variate  fitting  is,  in  fact,  accomplished  by  performing  a  series  of 
uni-variate  fittings.  The  initial  data  table  will  consist  of  the  following 
information: 


Xo 

X2 

^3 

^4 

^5 

Zo 

^00 

^01 

yo2 

yo3' 

yo4 

^05 

Zl 

^10 

^11 

^12 

^13 

^14 

^15 

Z2 

^20 

^21 

^22 

^23 

^24 

^25 

Z3 

^30 

^31 

^32 

^33 

^34 

735 

^4 

^40 

^41 

^42 

^43. 

744 

^45 

Zs 

^50 

^51 

^52 

^53- 

^54 

^55 

The  first  fitting  consists  of  the  representation  of  yQ  as  a 
function  of  X  at  Z  =  Zq.  The  result  is  an  equation  of  the  form: 


(8) 
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(9) 


■  ^00  ^01  ^  ^  ^02  ^  ^03  ^  ^04  ^  at  Z  -  Z^ 

This  process  is  repeated  at  Z^,  Z^,  Z^,  Z^,  and  Z^.  The  resulting  data 
can  then  be  written  as: 


xo 

xi 

X^ 

X^ 

X4 

^0 

^0 

^01 

A„o 

02 

^3 

^04 

Zl 

^10 

^1 

^12 

^13 

^14 

^2 

^20 

^21 

^22 

^23 

^4 

Z3 

^30 

^31 

^32 

^33 

^34 

Z4 

-^40 

A41 

^42 

^43 

^44 

Z5 

-^50 

^51 

^52 

^3 

A54 

The  resulting  coefficients  are  then  cross-fit  in  terms  of  powers  of  Z  in  a 
similar  fashion.  Thus,  the  first  column  of  line  (1)  results  in  an  equation 
of  the  form: 

*0  =  COpo  +  Djo  Z  +  °20  »  I>30  ^  X®  (U) 

This  process  is  repeated  for  each  of  the  other  columns.  The  result  is  an 
equation  of  the  form: 

y  "  Dqo  +  DqI  ^  *  ^02  ^  ^03  ^  ^04 

f  Djo  Z  f  Dii^Z  f  012^^2+  DjgX^Zf  0^4  Z 

+  D20  Z2  f  D2iXz  2  f  D22X.^z2  ^023  X^ZZ  f  D24  X^z2  (12) 

+  D30  Z3  f  D31  XZ3  r  D32  X2z3  +  D33  X^Z^  f  D34  X^ZS 

+  D40  Z4  f  D41  XZ4  4  D42  X2z4  f  D43  X3Z4  f  D44  X^Z^ 
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The  resulting  equation  (12)  must  then  be  unnormalized  to  the  form  of 
equation  (7).  The  process  Involves  a  series  of  linear  unnormalizing  steps, 
each  of  which  are  in  accordance  with  the  procedure  of  Table  3. 

C.  Three  or  More  Independent  Variables 

Those  elements  of  the  system  which  are  functions  of  three  or 
more  variables  can  be  represented  by  the  procedure  contained  In  this  section. 
The  result  will  be  an  empirical  equation  of  the  form: 

y  =  Boo  ^  ®01  ^1  Bo2  *2  *  ®03  *3  ■*■  •  •  •  ®0n  ^ 

f  Bi2  X1X2  Bi3  X1X3  f  . ..  +  Bln 

*■  hs  ^2*3  ^  •  +  B2n  ^2^^ 

f  Bi  xi^  f  B2  X2^  +  B3  X32  f  . . .  +  Bn  Xn^ 

where  xi  through  Xn  are  the  independent  variables  involved.  Note  that  no 
terms  of  order  three  or  greater  are  Included.  The  absence  of  these  higher 
order  terms  will  make  the  fitting  procedure  more  dependent  upon  the  user's 
ability  to  identify  the  proper  format  or  transformation  to  be  used  for  each 
variable. 

Although  no  specific  guidelines  can  be  offered  to  aid  In  the  identifica¬ 
tion  of  the  proper  transformations,  the  use  of  approximate  analyses  may  be 
extremely  helpful.  If  ineffective  transformations  are  used  for  any  of  the 
variables,  however,  this  factor  will  be  Indicated  by  the  results  of  the  fitting 


(13) 
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operation.  The  offender  could  then  be  replaced  by  a  more  suitable  variable 
transformation  and  a  new  fitting  obtained.  The  techniques  of  section  A  may 
be  helpful  in  determining  which  transforms  can  be  represented  with  sufficient 
accuracy  with  second-order  fittings  in  each  variable. 

After  the  variables  and  their  transforms  have  been  determined, 
an  experiment  plan  must  be  selected  to  match  the  number  of  variables  involved. 

The  number  of  coefficients  to  be  determined  in  developing  equation  (13)  is 
equal  to: 

Number  of  Coefficients  =  1  -f  2n  4-  — -  (14) 

2(n-2)! 

The  number  of  data  points  to  be  used  in  determining  the  coefficients  can  be 
found  from  the  equation: 

Number  of  Data  Points:  N  *  +  2n  4  k  (15) 

where  n  is  the  number  of  variables,  p  is  the  order  of  replication  (0  for  full 
replicate,  1  for  half,  2  for  quarter,  etc.),  and  k  is  the  number  of  center- 
points.  The  resulting  degrees  of  freedom  can  be  obtained  by  subtracting 
equation  (14)  from  equation  (15).  Table  4  contains  a  summary  of  these  data 
for  plans  involving  a  single  center-point. 

All  plans  must  have  at  least  one  degree  of  freedom  and  should, 
in  general,  have  5  to  10  degrees  of  freedom  in  order  to  maintain  a  balance 
between  the  number  of  points  required  and  the  number  of  terms  confounded  by 
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the  partial  replication  process.  Consequently,  the  following  plans  are 
recommended: 


equation: 


Full  Replicate  for  n  up  to  4 
Half  Replicate  for  n  of  5  and  6 
Quarter  Replicate  for  n  of  7 
Eighth.  Replicate  for  n  of  8 

The  orthogonal  factor  *  o<  -  must  then  be  calculated  from  the 


(16) 


Each  of  the 
-1  to  +  1. 


independent  variables  must  then  be  normalized  over  the  range  of 
This  is  accomplished  by  the  following  transform  equations: 
a  =  Wmin. 
b  =  (x)max. 


(17) 


(^)X  +  (j^) 


The  unnormalized  Independent  variables  must  then  be  calculated,  using  the 
last  of  equations  (17),  at: 

X  =  -o< 

=  -  1 

=  0  (18) 
=  1 

=  < 
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The  specific  experiment  plan  to  be  used  must  then  be  obtained  from  Table  5 
which  lists  the  specific  combinations  of  independent  variables  at  which  the 
dependent  variable  is  to  be  evaluated.  Table  5  also  lists  the  first  order 
interaction  terms  which  are  confounded  by  the  partial  replication  process. 

It  is  desirable  to  assign  independent  variable  numbers  so  that  only  marginally 
significant  effects  are  confounded.  Additional  experiment  plans  are  available 
from  Reference  (4) . 

The  dependent  variable  must  thdn  be  evaluated  at  all  levels 
indicated  by  the  experiment  plan.  The  results  of  line  (18)  will  assist  in 
converting  the  coded  values  of  the  independent  variables  in  the  experiment 
plan  to  actual  values. 

The  next  step  in  the  procedure  will  be  to  transfer  this  informa¬ 
tion  to  the  computer  input  sheet  which  will  contain  the  necessary  Instructions 
for  the  key  punch  operator.  Some  additional  Information  will  be  required  and 
this  will  be  Indicated  where  necessary.  A  sample  input  sheet  is  Illustrated 
in  Table  6. 

Section  A  contains  information  which  will  serve  to  identify  the 
specific  problem,  the  Investigator,  and  the  date  of  calculation. 

Section  B  contains  the  type  of  fractional  replication  selected 

from  Table  4. 

Section  C  contains  the  factorial  part  of  the  ejcperiment  plan 
selected.  It  may,  however,  be  omitted  v\4ien  a  full  replicate  is  used.  Line 
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3X1  contains  the  coded  values  of  the  first  independent  variable  -  the  f  1  or 
-  1  values.  If  more  than  22  values  are  required,  the  remainder  must  be 
continued  on  a  second  card  as  illustrated.  Similarly,  3X2  describes  the 
second  independent  variable,  3X3  the  third,'  etc.  Note  that  it  is  not 
necessary  to  indicate  the  cross  values  (  -  ^ )  or  center-points. 

Section  D  contains  the  values  of  the  dependent  variable 
corresponding  to  the  experiment  plan  selected  and  to  the  sequence  used  in 
Section  C.  All  factorial,  cross,  and  center-point  values  must  be  included. 
The  first  card  starts  with  the  identification  3Y1  and  contains  the  first  9 
values.  Each  succeeding  card  is  started  in  the  second  column  and  contains 
the  next  10  values  until  all  data  have  been  included. 

Section  E  contains  the  F-ratio  selected  from  Table  7  correspond¬ 
ing  to  the  number  of  degrees  of  freedom  associated  with  the  experiment  plan 
and  the  probability  that  the  deviations  observed  could  occur  by  chance  alonp. 
It  is  recommended  that  probabilities  of  1  to  5  %  should  normally  be  selected. 

Section  F  contains  the  boundaries  on  each  of  the  independent 
variables.  Line  3Z1  contains  the  uncoded  minimum  value  of  the  first  variable 
and  the  maximum  value.  Similarly,  line  3Z2  is  for  the  second,  3Z3  for  the 
third,  etc.  Note  that  these  values  must  correspond  to  the  -  1  and  f  1  coded 
values. 
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Section  G  contains  control  Information  as  follows: 


M  Is  the  number  of  Independent  variables  (10  maximum) 

N  Is  the  number  of  data  points 

NC  Is  the  number  of  center- points  (10  maximum) 

NR  Is  the  number  of  multiple  replications  (normally  1) 

NP  Is  the  order  of  fractional  replication  (0  for  full,  1  for 
half,  etc.) 

AN  Is  zero  If  analysis  of  variance  Is  desired,  otherwise  one. 
RSOP  Is  zero  If  response  surface  optimization  Is  not  desired, 
otherwise  one. 

HIGH  Is  zero  If  higher  order  Interactions  are  not  desired, 
otherwise  one. 

CONF  Is  zero  If  confidence  intervals  are  not  desired, 
otherwise  one. 

TPTS  Is  zero  if  the  unnormallzed  model  equation  Is  not  desired, 
otherwise  one. 

LOGE  Is  zero  If  a  model  based  upon  log  y  Is  not  desired, 
otherwise  one. 

ACT  should  be  Identical  with  TPTS 
It  Is  recommended  that  the  following  values  be  used  normally, 

NC  -  1  CONF  *  1 

NR  -  1  TPTS  «  1 


-14- 


AN 


LOGE  =  0 


=  0 

RSOP  =  0  ACT  -  1 

HIGH  =  0 

The  input  sheet  can  then  be  processed  and  the  data  run  on  the 
computer.  Results  should  be  checked  for  key  punch  errors  by  comparing  the 
print-out  against  the  input  sheet.  The  results  should  then  be  evaluated. to 
determine  if  the  fit  is  of  sufficient  accuracy  for  use  in  the  subsequent  system 
optimization.  The  following  data  will  be  printed: 

1.  The  uncoded  values  of  each  of  the  independent  variables 
at  the  -  df,  -  1,  0,  1,  and  levels.  These  data  should  correspond  with 
manual  calculations  using  equation  (18). 

2.  The  experiment  plan  which  shows  the  coded  values  of 
all  Independent  variables  at  each  of  the  data  points.  These  data  should 
correspond  with  the  plan  selected  and  with  the  data  of  section  C  of  Table  6. 

3.  The  dependent  variable  table  which  should  correspond 
to  section  D  of  Table  6. 

4.  The  comparison  of  observed  and  calculated  values  of 
the  dependent  variable  at  each  of  the  input  data  points.  In  addition  to  a 

direct  comparison,  several  measures  of  error  are  available.  The  most  significant 
of  these  is  the  term  labeled  PCT/RNGE  which  is  determined  from  the  formula: 


PCT/RNGE  = 


Calculated  y  -  Observed  y 
(Observed  y)  -(Observed  y)^t„^ 


(19) 


This  measure  indicates  the  efficiency  of  the  fitting  operation  since  a  maximum 
of  50  %  could  have  been  obtained  by  representing  the  model  as  a  constant. 
Although  some  judgment  is  peeded  in  determining  an  acceptable  error  criteria 
the  following  guidlines  may  be  used: 

0  to  2  %  is  an  excellent  fit 
2  %  to  10  %  may  be  an  acceptable  fit 
Above  10  %  is  generally  a  poor  fit 

If  a  poor  fit  is  obtained  with  correct  input  data ,  this  may  be  remedied  by 

using  an  alternate  transformation  for  one  or  more  of  the  independent  variables. 

+  , 

The  variables  which  should  be  transformed  are  those  which  are  at  their  -  (\ 
values  at  the  points  where  the  largest  errors  are  obtained.  No  specific  guide¬ 
lines  can  be  offered,  however,  to  assist  in  determining  the  proper  transforma¬ 
tion  to  be  used. 

5.  The  variance  of  the  estimate  and  coefficients.  The  term 
labeled  standard  error  of  the  estimate  is,  in  actuality,  the  variance  of  the 
estimate  and  is  a  measure  of  the  dispersion  of  the  individual  observations 
about  their  mean.  It  represents  another  measure  of  the  goodness  of  fit  of  the 
derived  analytical  model. 

6.  The  analysis  of  variance.  This  tabulation  contains  the 
F-ratlo  associated  with  each  of  the  derived  coefficients  of  the  model.  The 
highest  F-ratlo  corresponds  to  the  most  significant  term  in  the  model  and  pro¬ 
vides  a  means  for  ranking  the  individual  terms  in  prde^  of  their  significance. 
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All  terms  which  have  an  F-ratlo  less  than  the  value  in  section  E  of  Table  6 


can  normally  be  ignored. 

7.  The  coefficients  of  the  model  in  terms  of  the  coded 
values  of  the  Independent  variables.  The  number  in  parenthesis  indicates 
whether  or  not  the  coefficient  is  a  significant  one  within  the  prescribed 
value  of  F-ratio  -  1  indicates  the  term  is  significant  and  0  that  it  is  not. 

8.  A  repeat  of  line  4  using  only  the  significant  terms. 

If  the  measures  of  error  are  significantly  higher  than  those  obtained  in  line 
4,  the  selected  F-ratio  is  probably  too  high  and  some  other  value  should  be 
used. 

9.  The  coefficients  of  the  model  in  terms  of  the  uncoded 
values  of  the  Independent  variables  including  both  significant  and  insignifi¬ 
cant  terms.  Insignificant  terms  can  be  eliminated  by  referring  to  the  results 
of  line  7.  This  should  be  done  with  caution  since  the  unnormalizing  process 
will  result  in  contributions  to  the  linear  terms  from  the  quadratic  terms.  The 
results  of  this  last  step  may  be  utilized  directly  in  the  subsequent  system 
optimization  process. 
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V.  GLOSSARY 


Full  Factorial  Design  -  An  array  of  n  Independent  variables  (factors)  at  2,  3, 
or  more  coded  levels  of  each  variable  in  which  all  combinations  of  coded  levels 
and  factors  are  used  simultaneously  to  establish  the  characteristics  of  a  dependent 
variable.  The  number  of  treatments  or  combinations  required  is  2^,  3^,  etci,  All 
feasible  linear  effects  or  combinations  of  variables  can  be  determined  independently. 

Fractional  Factorial  Design  -  An  array  in  which  only  1/2 ,  1/4 ,  etc .  of  the 
full  factorial  design  is  used  to  establish  the  characteristics  of  the  dependent 
variable.  The  number  of  treatments  required  is  equal  to  2*^“^,  3”“P,  etc.  where 
p  is  the  order  of  replication  and  is  equal  to  0  for  a  full  replicate ,  1  for  a  half 
replicate/  2  for  a  quarter  replicate,  etc.  It  is  used  to  reduce  the  number  of 
treatments  required. 

Box-Wilson  Factorial  Design  -  A  modified  fractional  factorial  design  in  which 
the  2*^“P  treatments  are  increased  by  Jc  center  points  at  the  center  of  the  region  and 
2  n  cross  points  at  +  a  values  of  each  factor  successively.  The  additional  treat¬ 
ments  permit  the  estimation  of  quadratic  terms  for  each  variable. 

Confounding  -  The  process  by  which  the  linear  effects  are  grouped  in  order  to 
permit  estimation  with  a  fractional  factorial  design. 

Alias  Pairs  -  Specific  effects  which  are  confounded  by  the  partial  replication 
or  fractional  factorial  design.  Each  alias  pair  will  consist  of  two  effects  with 
a  half  replicate ,  three  with  a  quarter  replicate ,  etc . 
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Defining  Contrast  -  The  pattern  which  Identifies  the  specific  treatments 


retained  by  a  fractional  replicate  and  the  effects  which  are  confounded. 

Degrees  of  Freedom  -  The  difference  between  the  total  number  of  data  points 
or  treatments  used  and  the  number  of  coefficients  which  are  calculated. 

Variance  -  The  mean  value  of  the  sum  of  squares  of  the  deviations  of  a 
number  of  observations  about  their  mean. 

Standard  Deviation  -  The  square  root  of  the  variance . 

Objective  Function  -  The  variable  that  Is  to  be  optimized  subject  to  the 
constraints  provided. 

Constraints  -  Linear  or  non-llnear  boundaries  placed  upon  the  Independent 
variables  In  the  optimization  process . 
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Table  1 


Selection  of  the  Degree  of  Polynomial  Representation 

a.  Independent  Variable  Values 


Degree 

0 

1 

2 

3 

4 

^0 

0 

0 

0 

0 

0 

1.000 

.500 

.250 

.150 

.095 

^2 

1.000 

.750 

.500 

.345 

^3 

1.000 

.850 

.654 

1.000 

.904 

-5 

1.000 

Error  Coefficients 

Degree 

0 

1 

2 

3 

4 

Cq 

+  .5 

+  .25 

+  • 166 ... 

+  .125 

+  .10 

Cl 

-.5 

-.50 

-.333. . . 

-.250 

-.20 

^2 

+  .25 

+.333. . . 

+  .250 

+  .20 

C3 

- .166. . . 

-.250 

-.20 

C4 

+  .125 

+  .20 

-.10 

Table  2  (Continued) 
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Table  4 


Solution  of  Experiment  Plan 


Full 

Half 

Quarter 

Eighth 

lies 
— H 

Coefficients 

Points 

d.f. 

Points  d.f. 

Points  d.f. 

Points  d.f. 

■ 

3 

5 

2 

6 

9 

3 

10 

15 

5 

15 

25 

10 

21 

43 

22 

27  6 

28 

77 

49 

45  17 

29 

1 

36 

79  43 

47 

11 

45  , 

81 

36 

49  4 

55 

147 

92 

83  28 

66 

149  83 

Table  5 


DESIGN  MATRIX  PLAN 


1 .  Number  of  independent  variables  -  n  «  3 

2.  Full  Replication  (p  ■  0) 

3.  Number  of  data  points  -  N  ■  (2’'^  +  2  n  +  K)  where;  K  ■  number  of  center  points 

4 .  Alpha  value  a  is  obtained  from: 

a  -  1/2  (  2"-‘>  (2”-''  +  2n  +  K)  -  2"-'>P 


N  Number  of 
Observations 

’‘i 

==2 

Factor  Level 

r  1 

1 

1 

1 

2 

1 

1 

-1 

3 

1 

-1 

1 

Full  Factorial 

4 

1 

-1 

-1 

Portion 

< 

5 

-1 

1 

1 

6 

-1 

1 

-1 

7 

-1 

-1 

1 

L  8 

-1 

-1 

-1 

f  ® 

a 

0 

0 

10 

-a 

0 

0 

6  Alpha  Values 

11 

0 

a 

0 

< 

12 

0 

-a 

0 

13 

0 

0 

a 

.14 

0 

0 

-a 

Center  Point 

{15 

0 

0 

0 
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DESIGN  MATRIX  PLAN 

1 .  Number  of  Independent  variables  -  n  ■  4 

2 .  Full  Replication  (p  »  0) 

3.  Number  of  data  points  -  N=(2  ^  +  2n-»-K)  where:  K  ■  number  of  center  points 

4 .  Alpha  value  -  a  is  obtained  from: 

a  »  1/2  (  2^“^  (2'^'*’  +  2  n  +  K)  -2""^) 


N  Number  of 
Observations 

^1 

Factor  Level 
^3  ^4 

» 

1 

1 

1 

1 

1 

2 

1 

1 

1 

-1 

3 

1 

1 

-1 

1 

4 

1 

1 

-1 

-1 

5 

1 

-1 

1 

1 

6 

1 

-1 

1 

-1 

7 

1 

-1 

-1 

1 

Full  Factorial 

1 

8 

1 

^1 

-1 

-1 

Portion 

9 

-1 

1 

1 

1 

10 

-1 

1 

1 

-1 

11 

-1 

1 

-1 

1 

12 

-1 

1 

-1 

-1 

13 

-1 

-1 

1 

1 

14 

-1 

-1 

1 

-1 

15 

-1 

-1 

-1 

1 

1  16 

-1 

-1 

-1 

-1 

17 

a 

0 

0 

0 

18 

-a 

0 

0 

0 

19 

0 

a 

0 

0 

8  Alpha  Points 

20 

0 

-a 

0 

0 

21 

0 

0 

a 

0 

22 

0 

0 

-a 

0 

23 

0 

0 

0 

a 

24 

0 

0 

0 

-a 

Center  Point 

1 

C25 

0 

0 

0 

0 
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Table  5  -  (Continued) 


DESIGN  MATRIX  PLAN 

1 .  Number  of  Independent  variables  -  n  =  5 

2 .  Half  Replication  (p  ■  1) 

3 .  Number  of  data  points  -  N  =  (2  *^  +  2n  +  K)  where:  K  =  number  of  center  points 

4 .  Alpha  value  -  a  is  obtained  from: 

a  =  1/2  (  2""^  (2"'^  +  2  n  +  K  -  2”'**) 


N  Number  of  ;■ 
Observations 

*2 

Factor  Level 
==3  ==4  ’‘S 

1 

1 

1 

1 

1 

-1 

2 

1 

1 

1 

-1 

1 

3 

1 

1 

-1 

1 

I 

4 

1 

1 

-1 

-1 

-1 

5 

1 

-1 

1 

1 

1 

6 

1 

-1 

1 

-1 

-1 

7 

1 

-1 

-1 

1 

-1 

Half  Factorial 

8 

1 

-1 

-1 

-1 

1 

Portion 

9 

-1 

1 

1 

1 

1 

10 

-1 

1 

1 

-1 

-1 

11 

-1 

1 

-1 

1 

-1 

12 

-1 

1 

-1 

-1 

1 

13 

-1 

-1 

1 

1 

-1 

14 

-1 

-1 

1 

-1 

1 

15 

-1 

-1 

“1 

1 

1 

ll6 

-1 

-1 

-1 

-1 

-1 

f  17 

a 

0 

0 

0 

0 

18 

-a 

0 

0 

0 

0 

19 

0 

a 

0 

0 

0 

20 

0 

-a 

0 

0 

0 

10  Alpha  Values 

1 

21 

0 

0 

a 

0 

0 

22 

0 

0 

-a 

0 

0 

23 

0 

0 

0 

a 

0 

24 

0 

0 

0 

-01 

0 

25 

0 

0 

0 

0 

a 

.  26 

0 

0 

0 

0 

-a 

Center  Point 

CM 

0 

0 

0 

0 

0 

Defining  Contrast: 
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Table  5  -  (Continued) 


DESIGN  MATRIX  PLAN 

1.  Number  of  variables  -  n  «  6 

2 .  Half  Replication  (p  ■  1) 

3 .  Number  of  data  points  -Ns  (2^^  ^  +2  n  -i-  K)  where:  K  «  number  of  center  points 

4 .  Alpha  value  -  a  is  obtained  from: 

a  -  1/2  (  (2""*’  +  2  n  +  K)  -  2”"*^) 


N  Number  of 
Observations 

=^1 

=‘2 

Factor  Level 
*3  *4  *5 

Jis. 

1 

-1 

-1 

-1 

-1 

-1 

-1 

2 

1 

1 

-1 

-1 

1 

1 

3 

1 

-1 

1 

1 

1 

-1 

4 

-1 

1 

1 

1 

-1 

1 

5 

1 

1 

-1 

-1 

-1 

-1 

6 

-1 

-1 

-1 

-1 

1 

1 

7 

1 

-1 

1 

1 

-1 

1 

8 

■rl 

1 

1 

1 

1 

-1 

9 

1 

-1 

1 

-1 

-1 

-1 

10 

-1 

-1 

-1 

1 

1 

-1 

11 

1 

1 

-1 

1 

-1 

1 

12 

•  -1 

1 

1 

-1 

1 

1 

13 

-1 

1 

1 

-1 

-1 

-1 

14 

-1 

-1 

-1 

1 

-1 

1 

15 

1 

-1 

1 

-1 

1 

1 

Half  Factorieil 

< 

16 

1 

1 

-1 

1 

1 

-1 

Portion 

17 

1 

-1 

-1 

-1 

1 

-1 

18 

-1 

1 

-1 

-1 

-1 

1 

19 

-1 

-1 

1 

1 

-1 

-1 

■ 

20 

1 

1 

1 

1 

1 

1 

21 

1 

-1 

-1 

-1 

-1 

1 

22 

-1 

1 

-1 

-1 

1 

-1 

23 

1 

1 

1 

1 

-1 

-1 

24 

-1 

-1 

1 

1 

1 

1 

25 

1 

-1 

-1 

1 

-1 

-1 

26 

-1 

-1 

1 

-1 

1 

-1 

27 

1 

1 

1 

-1 

-1 

1 

28 

-1 

1 

-1 

1 

1 

1 

29 

-1 

1 

-1 

1 

-1 

-1 

30 

-1 

-1 

1 

-1 

-1 

1 

31 

1 

1 

1 

>1 

1 

-1 

(32 

1 

-1 

-1 

1 

1 

1 

Continued  on  next  page 
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Table  5  -  (Continued) 


DESIGN  MATRIX  PLAN  -  Continued 

1 .  Number  of  variables  -  n  =  6 

2 .  Half  Replication  (p  »  1) 


N  Number  of 
Observations 

*2 

Factor  Level 
^3  ^4  ^5 

33 

a 

0 

0 

0 

0 

0 

34 

-a 

0 

0 

0 

0 

0 

35 

0 

a 

0 

0 

0 

0 

36 

0 

-a 

0 

0 

0  • 

0 

37 

0 

0 

a 

0 

0 

0 

38 

0 

0 

~<i 

0 

0 

0 

12  Alpha  Values 

39 

0 

0 

0 

a 

0 

0 

40 

0 

0 

0 

-a 

0 

0 

41 

0 

0 

0 

0 

a 

0 

42 

0 

0 

0 

0 

-a 

0 

43 

0 

0 

0 

0 

0 

a 

44 

0 

0 

0 

0 

0 

-a 

Center  Point 

(45 

0 

0 

0 

0 

0 

0 

Defining  Contrast:  X^X2X2X^XgXg 
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Table  5  -  (Continued) 


DESIGN  MATRIX  PLAN 

1 .  Number  of  variables  -  n  «  7 

2.  Quarter  Replication  (p  =  2) 

3.  Number  of  data  points  '-Nae(2  +2n  +  K)  where:  K  ■  number  of  center  points 


Continued  on  next  page 
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Table  5  -  (Continued) 


DESIGN  MATRIX  PLAN  -  Continued 

1 .  Number  of  variables  -  n  ■  7 

2 .  Quarter  Replication  (p  =  2) 


N  Number  of 
Observations 

*2 

Factor  Level 
*3  *4  ’‘S 

==6 

*7 

'33 

0 

0 

0 

0 

0 

0 

34 

'  -a 

0 

0 

0 

0 

0 

0 

35 

0 

a 

0 

0 

0 

0 

0 

36 

0 

-a 

0 

0 

0 

0 

0 

37 

0 

0 

a 

0 

0 

0 

0 

38  ' 

0 

0 

-a 

0 

0 

0 

0 

14  Alpha  Values  * 

39 

0 

0 

0 

a 

0 

0 

0 

40 

0 

0 

0 

-a 

0 

0 

0 

41 

0 

0 

0 

0 

a 

0 

0 

42 

0 

0 

0 

0 

-a 

0 

0 

43 

0 

0 

0 

0 

0 

a 

0 

44 

0 

0 

0 

0 

0 

-a 

0 

45 

0 

0 

0 

0 

0 

0 

a 

146 

0 

0 

0 

0 

0 

0 

-a 

Center  Point  ‘ 

C47 

0 

0 

0 

0 

0 

0 

0 

Defining  Contrasts:  X,X„X,X.X_ 

1  2  3  4  5 

X.X-X-X 

4  5  6  7 

V6=V7'  Vs  ^ 

X5X,,  : 

K.  X,  =  X. 
4  7  5 

X. 

0 

- 

Alias 

Pairs 
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Table  6 


CQ 

<n 

o 

>-A 


<U 

CO 

(0 

M 

B 

0> 

•H 

(0 

I 


B 

O 

o 

V 

CO 


B 

g 

3 

O 

U 


S' 

A 

u 

B 

3 

a 

a 

<0 


w' 

ii 


(0 

o 

0)  Q 

C  CT)| 
£>  to 


CM 

B 

6 

5 

o 

O 

B 

2 

B 

3 

a 


3i 


®l 

B 

O 

a 

o 


B 

O 

u 

a 

a 

<s 

3 

■fO 

P 


9 


m 


I  I  I  I  I  I 


I  . 

I  I  I  I  I 
.  H  I  H  M 

CM|  |co|  IM*! 

XI  1X1  1X1 

c*>|  |co|  |co| 


H 

ml 

XI 

col 


I 

H 

<3| 

XI 

co| 
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1 

i 


1 


1 


NC  NR  NP  AN  RSOP.  HIGHCQl^T  TPTS^LOGC'^T 


Degrees  of 
Freedom 


TABLE  7 


F-Ratlo 


161.4 

18.51 

10.13 

7.709 

6.608 


5.987 

5.591 

5.318 

5.117 

4.965 


4.844 
4.747 
4.667 
.600 
.543 


Probability 

.01 


4052 

98.50 

34.12 

21.20 

16.26 


3.74 

2.25 

1.26 
10.56 
10.04 


9.646 

9-330 

9.074 

8.862 

8.683 


405,284 

998.5 

167.5 
74.14 
47.04 


35.51 

29.22 

25.42 

22.86 

21.04 


9.69 

8.64 

7.81 

7.14 

6.59 


4.325 

4.301 

4.279 

4.260 

4.242 


8.017 

7.945 

7.881 

7.823 

7.770 


14.59 

14.38 

14.19 

14.03 

13.88 
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